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Abstract. We review some recent results concerning Gibbs measures for nonlinear 
Schrodinger equations (NLS), with implications for the theory of the NLS, including 
stability and typicality of solitary wave structures. In particular, we discuss the 
Gibbs measures of the discrete NLS in three dimensions, where there is a striking 
phase transition to soliton-like behavior. 



1. Introduction 

The cubic nonlinear Schrodinger equation (NLS) arises naturally in physics and 
engineering, for instance, as a model of nonlinear optics, Langmuir waves in plasmas, 
and Bose-Einstein condensation (BEC) J3U [5IJ HH9 HUB EH ■ It turns out that the NLS 
can be rigorously derived from the physics of interacting particles in BEC via a scaling 
limit and statistical mechanics [201 ESS • It a ls° turns out that the probabilistic tools 
of statistical mechanics are key to understanding the dynamics of nonlinear Schrodinger 
systems, especially why solitons are so common in certain regimes. (They are so highly 
probable that they are the only structures observed macroscopically.) 

Consider a wavefunction u, that is, u mapping a time variable f £ 1 and space 
variable x € M. d into the complex plane C, which satisfies a d-dimcnsional nonlinear 
Schrodinger equation (NLS): 

(1.1) id t u = -Ait + k\u\ p ~\ 

where A is the Laplacian operator in M d , p is the nonlinearity parameter, and k is a 
coupling parameter which is either +1 (repulsive/defocusing interaction) or —1 (attrac- 
tive/focusing). As a matter of notation, we will usually write a continuous wavefunction 
in the form u(t,x) and a discrete wavefunction in the form fk(t) = f(t, k) where k is 
the discrete spatial variable. 

We will be interested in chiefly the cubic (p = 3) and focusing (« = —1) NLS, and 
especially its discrete counterpart via the finite difference approximation, where the 
spatial variable k is in a discrete set. Also of interest is a modification of ( |1.1[ ) with a 
fractional power of the Laplacian as a non-local dispersive operator, an equation arising 
naturally in biophysics models of macro- molecules like DNA (211 [25j |46j . 

The plan of this article is largely chronological, and background material for the 
continuous and discrete NLS can be found in [H [21 [55] ■ In section [2] we will discuss 
solitons and long-time results for the discrete NLS based on analogues of dispersive 
PDE methods. Then in section [3] there are similar but more nascent results for discrete 
NLS-type systems with long-range interactions, again adapted from PDE methods. We 
introduce the Gibbs measures for continuous NLS in section [4 and discuss applications 
to questions of well-posedness for low-regularity initial data (4.11), applications to the 
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typicality of solitons in the absence of wave collapse (4.2 ), and a controversial conjecture 



for the Gibbs measures of the NLS (4.3). There is a way to get around the problem 
of nonexistence of Gibbs measures in the presence of wave collapse, so in section [5] we 
discuss new results for the Gibbs measures of the 3D discrete NLS, where there is a 
striking phase transition to soliton-like behavior. 

2. Solitons and breathers for the discrete NLS 

Consider a lattice V which is 7L d , for the present, with lattice spacing h, i.e., a 
positive real number denoting the "distance" between two neighboring vertices in V 

The discrete NLS on V is a family of coupled ODEs governing the evolution of 
discrete wavefunctions /& = fk(t), i.e., maps from k 6 V and f 6 K into C satisfying: 

(2.1) *^/ fc = -A/ fc + K |/ fc r 1 / fc , keV. 

at 

Here the discrete nearest-neighbor Laplacian A is defined by 

:=i £(£"/*)• 

j~k 

Note that the sum is over all vertices j that are neighbors of k, and that division by h 2 
means that the discrete Laplacian converges to the continuous Laplacian as h goes to 
zero. 

Discrete NLS equations can exhibit localized modes or discrete breathers, which are 
periodic-in-time and spatially localized solutions. Breathers have a long history in the 
literature, and Flach, Kladko, and MacKay made a physics conjecture that the energies 
of breathers must be bounded away from zero in supercritical lattice models (e.g., 
discrete NLS in 3D) [23]. Weinstein proved one of the conjectures in the supercritical 
case, that ground state standing waves exist if and only if the total power is larger 
than some strictly positive threshold [50], called the excitation threshold, which is 
characterized by variational methods. 

In order to make this precise, we define the total power or particle number of a 
wavefunction / solving the discrete NLS ( |2.1[ ) as: 

|2 



N(f) := ||/||| = £ \f, 



fcl 
kev 

Also the Hamiltonian energy is well-known to be conserved in time: 
H(f) : = _(A/J) - -A- £ l/^r 1 . 

Here the inner product (•,•) is the usual one on l 2 (Z d ). The ansatz for breathers, or 
standing waves, is that they are of the form for oj £ K: 

/*(*) = e~ iut g x , k&V, f k (t)el 2 (V). 

Then g solves the discrete Euler-Lagrange equation 

(2.2) wg x = -Ag x - \g x \ p ~ 1 g x . 



And the solution of (2.2), called the ground state, can be constructed as a minimizer 
of a variational problem: 

I v = mf{H(f) : N(f) = u}. 
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Then for I v < 0, the minimum is attained by g (which is not necessarily unique), with 
some ui = oj(v) such that both the constraint N(g) = v and the Euler-Lagrange system 



(2.2) hold. Then the main result of [60] is the following. 



Theorem 2.1 (Weinstein [60] L For supercritical nonlinearities p > g + 1 (for example, 
cubic nonlinearity in three dimensions), there exists an I 2 -excitation threshold v c = 
v c (d,p) > for the ground state of the discrete NLS, meaning that: 

l v < if and only if v > v c . 

On the other hand, for subcritical 1 < p < | + 1 (cubic in one dimension), there is 
no excitation threshold, i.e., I v < for all v > 0. 

As a preview of the statistical mechanics perspective we will take later, Weinstein's 
result means that there exist localized modes at zero temperature, corresponding to 
the ground state. (In the context of statistical mechanics, temperature is a scaling 
parameter that determines how likely states of various energies are.) Using statistical 
mechanics, we will be able to add that such localized modes exist for a whole range of 
temperatures and moreover are typical in a probabilistic sense. Moreover, there will 
be a phase transition as the mass and the temperature are varied: For low values of 
mass or inverse temperature, the states that are highly probable are small states with 
negligible energies. But then there's a sudden jump at a critical threshold of mass or 
inverse temperature to highly probable localized states with negative energies, in fact 
bounded away from zero. 

Related work of Malomed and Weinstein used this same variational approach to an- 
alyze the dynamics of solitons in the ID discrete NLS: For nonlinearities up to about 
cubic, there is a unique soliton of large width and no excitation threshold. For nonlin- 
earities in a medium range between cubic and quintic, there is no threshold and there 
are three solitons, one unstable and two stable (in sharp contrast to the situation for the 
continuous NLS). And for nonlinearities above quintic there is an excitation threshold 
for the solitons, one of which is stable and highly localized, and the other unstable and 
wide [33] . 

Weinstein also made a conjecture in [60] about the subthreshold behavior, that if 



N{f) < v c for the solution / of the discrete NLS (2.1 1, then / disperses to zero, in the 
sense that: 

\\f{t)\\l«{v) ^* 0, for any q e (2, 00]. 

This conjecture has been addressed by similar adaptations of PDE methods to the dis- 
crete NLS by Stefanov and Kevrekidis, who proved the following decay and Strichartz- 
type estimates for discrete Schrodinger (and Klein-Gordon) equations. 

Theorem 2.2 (Stefanov-Kevrekidis 57]). For the free discrete Schrodinger equation 
i4ffk = —Afk with initial condition /(0) = jfc(0), we have the decay estimate for the 
wavefunction f(t) = fk(t): 

||/(t)||j- <C(l + |t|)- d / 3 ||/(0)|| ;1 . 



And for the nonlinear discrete Schrodinger equation (2.1), with initial data f(0) in the 
sequence space I 2 , we have the Strichartz-type estimates for any admissible pairs (q,r) 
and (q, f) (admissible means 2 < (q, r) < 00, | + ^ < | , and (q, r, d) ^ (2, 00, 3) ), that 

ll/Wlk^<c(||/(o)|| ;2 + |||/| 2 /ll^^), 
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where q' is the dual exponent of q, and the mixed Lebesgue spaces L q t V x are defined by: 
ll/*(<)||w:= ( r(j2\^\ r ) dt 



kev 



Moreover, the constant in the bound can be chosen uniformly in time for any compact 
interval [0, T]. 

The decay here is not as fast as for the continuous Schrodinger decay estimates 
(which is of order t~ d / 2 ), but it is sharp in the discrete setting. One consequence of this 
approach is the following result about decay of small solutions, addressing Weinstein's 
conjecture and showing that standing wave solutions are impossible in the subthreshold 
regime. 

Corollary 2.1 (Stefanov-Kevrekidis [57]). If p > § + 1, then there exists a positive 
e depending on dimension only, such that for small enough initial data, ||/(0)||;2 < e, 
a unique solution exists for all time and vanishes in V in the limit as t —¥ oo, for all 
r > 2. In fact: 

\\f\\ L nr<Ce. 

Moreover, under certain conditions (p > 4 + 1 and d < 2), small solutions decay like 
the solution to the free Schrodinger equation above: for all q G [2, (8 — 2d)/(d + 1)], 

\\f(t)\\i« <^- d(g - 2)/(39) ll/(o)IU<. 

This is a scattering result, where the solution of the nonlinear equation is asymptotic 
to (scatters to) the solution of the linear (free) Schrodinger equation as time goes to 
infinity |58) . But this uniform smallness of the solution for small initial data happens 
only when the nonlinearity is more than just supercritical (% + 1 rather than § + 1). 

There are two scaling limits of the discrete Schrodinger equation that are of interest. 
One is the continuum limit alluded to earlier, of sending the lattice spacing h to zero. 
In this limit, at least in one and two dimensions, one obtains the continuous nonlinear 



Schrodinger equation ( 1.1 1. Another limit is the anti-integrable or anti- continuum limit, 
resulting in an infinite system of decoupled ODEs by sending h to infinity, for instance, 
or equivalently sending the total power v to infinity. In this case, ground states become 
concentrated on the lattice at a single site [5U] , and their spectral properties and stability 
issues have been studied extensively [36 , 33, 49, 50, 51j. 

(An interesting thing will happen for the statistical mechanics of the continuum limit 
in three dimensions, however, because the nonlinearity is supercritical there and thus 
dominates the kinetic energy, making the coupling effectively negligible and producing 
behavior similar to the anti-continuum limit in that states typically localize at a single 
site.) 

For instance, Kevrekidis, Pelinovsky, and Stefanov |33j build on the previous disper- 
sive decay estimates of [FT] , in order to prove asymptotic stability of the ground states 
under some assumptions on the spectrum of the linearization of the DNLS equation 
with external potential in ID: 

(2.3) *^A = (-A + Vfc)/ fc + K |/ fc r 1 /fc, keV, 

Because of the slower rate of decay in the l°° norm of solutions to the discrete Schrodinger 
equations \2.2\ the critical power of the nonlinearity is p = 7 in place of the continuous 
Schrodinger equation's critical power p — 5, and [S3] does not handle the cubic case 
p = 3. 
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Define the spaces l] and I 2 by the norms: 

H/11,1 :^^(l + fc 2 ) s / 2 |M, 
fcez 

\\f\k ■= (E( 1 + fc2 ) s i^i 2 

Vfeez 

Decompose the solution into a family of solitons and a radiation part: 

<t) = e-^Ht)) + z(t)]. 

Theorem 2.3 (Kevrekidis-Pelinovsky-Stefanov [33]). Assume that k = 1, that p > 4, 
that V £ l\ a for some a > |, t/iai no solution ipo of Htpo = exists in l 2 _ s for 
s £ (5, |], and £/ia£ := —A + V has exactly one negative simple eigenvalue usq < 
with eigenvector t/jQ £ I 2 and no eigenvalues greater than 4. Also assume that there is 
no initial phase shift 8(0) — 0, that the initial soliton location w(0) is close to a fixed 
luq, and that the initial data u(0) is close to the soliton <fi(uj(0)) in the space I 2 . 

Then in the limit as t goes to infinity, uj(t) converges, and the radiation part z(t) 
vanishes in l°° . 

This is a nice asymptotic stability result, but what's more remarkable is that some- 
thing similar happens in a continuum limit of high-dimensional (d > 3) finite lattices 
called the thermodynamic limit, sending h to zero and simultaneously taking the total 
mass (power divided by number of lattice sites) to be fixed (i.e., sending power to in- 
finity in particular way). In the thermodynamic limit, one can find a critical threshold 
above which one obtains states concentrated states at a single site, at least with high 
probability. 

3. SOLITONS FOR DISCRETE NLS WITH LONG-RANGE INTERACTIONS 

In pure physics applications like optics and Bose-Einstein condensation, the NLS is 
the model, but in the context of biophysics where particle interactions can be long- 
range, the natural model turns out to be a Schrodinger equation with fractional-order 
Laplacian. An example is modeling electron transport in DNA j24j [25J 06] . On DNA, 
and other macromolecules with complicated helical structures, electrons or excitations 
can propagate along the strand in a very non-local way. 

Mathematically, these phenomena are just beginning to be explored. The usual ques- 
tions about well-posedness and solitons are interesting for the fractional NLS, and in 
fact, there is a corresponding discrete system whose solitons have been studied numeri- 
cally and whose continuum limit is the predicted fractional NLS. Work on the fractional 
PDEs has been done by jT3J [27] , on fractional NLS specifically by [3J [22] , and on the 
related Benjamin-Ono equation by [321 H5] , 

The discrete set-up is on a ID lattice KZ with mesh size h > 0, and some fixed ho 
such that h < ho ^ 1. Writing Xk — hk with k £ Z, we consider a family of discrete 
wavefunctions f£ — f h (t, Xf.) : K x KL — > C that satisfy the following discrete NLS-type 
equation: 

(3-D (*) = h £ §^§M + "itfwrvft*) 

Here the sum is over all vertices j distinct from fixed k, and < s < 00 is a fixed 
parameter governing the decay of the long-range interaction. These wavefunctions 
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have conserved quantities similar to the traditional discrete NLS for energy: 

TT(fh\_ fl2 l/fc~/j| 2 Kh V-^|ffe|P+l 



and for mass: 



M(f h ) = hJ2\fk\ 



k 

Equivalent to this last, power is conserved: 



N{f h ) = Y,\ti\ 



As a model of DNA, the cubic nonlinearity describes a self-interaction for a base 
pair of the strand with itself, and the summation term models interactions between 
base pairs decaying like an inverse power of the distance along the strand [46 . The 
complex twisted structure of DNA in three dimensions is what allows distant base pairs 
to interact, even at a long distance. 

Under reasonable conditions, solutions of this "fractional" discrete Schrodinger equa- 



tion (3.1 1 exist for all time and are unique by standard fixed-point methods [31 l2"2l |3~4"] . 
and solitary wave solutions of this system of highly coupled ODEs have been derived 
[2"1] and extended to more general nonlinearities [3]. 



Explicit solitary wave solutions for this discrete fractional NLS (3.1 1 have been stud- 
ied in the biophysics literature. For instance, Gaididei et al. [24] make the ansatz 
f h (xk) = e lXt <l> h (xk), w ith localized state: 



, h, % h 1 iVsinha , ., ... 

(3-2) ^k) = ^ cQsHa{2S _ 1)) e* P (-a\k-5\), 

the parameter a representing the inverse "width" of the soliton and i5 € [0, 1) de- 
termining its location. Numerically studying the associated variational problem, and 
minimizing H under the power constraint N(f h ) — constant, they find two kinds of 
localized states, one called "on-site" with 6 = and one "inter-site" with 8 = 1/2, 
depending on the value of the inverse power parameter s. 

For on-site states with cubic nonlinearity, there appears to be a critical value s c , esti- 
mated to be near 1, above which the dispersive interaction decays fast enough that the 
behavior is qualitatively like that of the nearest-neighbor approximation, or traditional 
discrete NLS. And below the critical value s c , they find three bistable localized states, 
one low-frequency wide state and two high-frequency narrow states [24j . On the other 
hand, inter-site states with cubic nonlinearity are not so well understood, and they find 
a lower critical value s' c ~ .55 and a smaller regime s € (2, s' c ) in which there are two 
stable inter-site states. In both cases, the dependence of the power on the frequency of 
the localized states is studied. 

Gaididei et al. also present the continuum version of these nonlocal models, with 
u = u(t, x) solving the fractional NLS of the form: 

(3.3) id t u = c(-A) a u±\u\ p - 1 u 

with u : K x K — > C, a constant c depending only on s, and a depending on s ap- 
propriately. In the biophysics literature, this limit h — > is routinely taken without 
mathematically rigorous arguments, but rigor is forthcoming with the result as follows 

M. 
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We consider the cubic interaction p — 3 for simplicity, but this can be generalized to 



other nonlinearities p, and we consider the inverse power law decay in (3.1), but this 
can be generalized to other interaction kernels in [34J that fall into asymptotic classes 
with behavior substantively similar to the kernel in with parameter s. 

The result is that if we are given a family of initial data for the discrete evolution 



problems (3.1), < h < ho, that approximate a sufficiently smooth initial datum for 



the fractional NLS (3.3), then the solutions of the discrete equations converge in an 
appropriate weak sense to the solution of the fractional NLS, with the fractional power 
of the Laplacian: 



a(s) 



s, for | < s < 1, 
1, for s ^ 1. 



In other words, the critical value is exactly s c = 1 (not too far from the numerical re- 
sults in [21]) below which the long-range interactions produce a nonlocal fractional NLS 
in the continuum limit, with Laplacian of order a = s. Above the critical threshold, the 
interaction strength decays so quickly that only local effects survive in the continuum 
limit, which is exactly the "classical" local NLS, a = 1. And at the threshold, we get 
the classical NLS in the continuum limit, with a logarithmic factor appearing in some 
scaling constants. 

For this continuum limit result, we develop fractional discrete Sobolev-type inequal- 
ities in the discrete setting, uniform embedding, interpolation, and a priori estimates, 
and finally for the limit itself, suitable weak compactness arguments. It would be nice 
also to obtain rigorous results about the solitary wave states from the biophysics liter- 
ature and their stability for both the discrete and continuous long-range systems, and 
these long-range NLS systems should be amenable to the statistical mechanics approach 
discussed next. 

4. Statistical mechanics and Gibbs measures for the continuous NLS 

Adapting tools from the theory of PDEs like variational methods and Strichartz 
estimates has been useful for understanding the discrete NLS and its localized soliton 
states, but it is nice to get beyond standard existence results to generic behavior. By 
setting up the Gibbs measure, or statistical ensemble, which is a probability measure 
on the space of all solutions, we can tell which behaviors are typical in a probabilistic 
sense. 

Gibbs measures, or invariant measures, are based on normalizing the Boltzmann 
distribution, so that the probability of finding the system in state u is: 

1 r -0H(u) 



Here f3 is the inverse temperature, sometimes replaced by (fc^T) -1 , which makes tem- 
perature T explicit using the Boltzmann constant kg. And Z(f3) is the normalizing 
constant, also known as the partition function, which encodes statistical properties of 



the system (in particular the free energy as we'll see in section 5.1). For dynamical 
systems such as the NLS, the energy H is conserved by the dynamics, so that we can 
invoke the Liouville theorem to see that the measure is invariant. 

The intuition behind the Boltzmann distribution can be understood by thinking 
about the extreme cases: In the high temperature case, (3 -C 1, the Gibbs measure is 
essentially uniform, that is, the system could statistically be in almost any state. On 
the other hand, near zero temperature, /? > 1, the Gibbs measure is nearly vanishing 
for any state but the lowest-energy one, and thus the system has a strong preference for 
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the ground state(s) . The statistical ensemble associated to the Gibbs measure is called 
the canonical ensemble, and the states that the system is likely to be in, according to 
the Gibbs measure, are called canonical macrostates, so called because those are the 
states that are observed macroscopically. 

There are variations on the standard Gibbs measure and its canonical ensemble. For 
instance, the microcanonical ensemble arises from conditioning the Gibbs measure on 
the energy H being a constant E (really being within e of E and sending e to zero 
appropriately). A natural question to ask is whether the microcanonical ensemble is 
equivalent to the canonical ensemble, meaning that they share the same macrostates. 
The answer is yes and no: if certain conditions are satisfied, they are equivalent, and 
if not, the microcanonical ensemble has a richer set of macrostates. Another ensemble 
is the petit canonical ensemble, which results from conditioning on mass: this is also 
called a mixed ensemble because it is canonical in energy and microcanonical in mass. 

An advantage of Gibbs measures is that they ignore transient states and focus instead 
on the long-time probabilistic behavior of the system, and this approach was initiated 
in the work of Lebowitz, Rose, and Speer [39] for the ID focusing NLS on the circle T: 

(4.1) id t u = -d xx u - |u| p ~V 

The canonical ensemble isn't the right thing to look at because it is not normalizable, so 
they conjectured that, with a mass cutoff, one can construct a formal invariant measure 
from the Hamiltonian energy: 

(4.2) H{u) := \ ( L \u'(x)\ 2 dx -L. t \ u {x)\^dx. 

2 Jo P + 1 J a 

The reason for the mass cutoff is that if one forms the Gibbs measure with density 
e -f3H(u)^ t nen ^ j s not normalizable because Z(f3) is infinite. This can be fixed by 
making B the allowed mass: the partition function Z is then a (finite) function of both 
the thermodynamic /3 and the allowed mass B. Thus we form the following formal 
invariant measure: 

(Z^B^e-^h^^yHduix). 

x 

Here the product needs to be interpreted properly, highlighting one of the difficulties 
in constructing Gibbs measures for infinite-dimensional systems: a finite-dimensional 
approximation must be done, for instance, by spectral truncation or by spatial dis- 
cretization. 

4.1. Application of Gibbs measures to well-posedness. Bourgain continued the 
investigation by using spectral truncation to avoid the harmonic analysis on cyclic 
groups that would have been necessary for spatial discretization [4]. He introduced 
a discretization in Fourier space by truncating the Fourier sum at the n-th Fourier 
coefficient: 

u n := P n u := J2 u(k)e 2mkx . 

|fe|<n 

Then the finite-dimensional Hamiltonian system is: 

(4.3) id t u n = -d xx u n - P n (|uT _1 , 
with the Hamiltonian energy functional: 

p+i 



life, p + lJ ^ 



E 2-nikx 
a k e 

|fe|<n 



dx, 
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where a,k is the fc-th Fourier coefficient: au '■= u(k). The Hamiltonian is conserved, as 
is the truncated mass: 



Next consider a ball in this finite-dimensional phase space: 



Sir. 



\k\^n 



\a k \ 2 < i? 2 



and the natural Gibbs measure on tt n .B at fixed inverse temperature ft — 1 is then: 

p+i 



dfl r , 



1 



exp 



P + 1 A 



E„ Jl-xikx 
afce 

|fc|<n 



cte >p„ dag, 



where /0„ := Z n 1 exp{— 27r 2 X)o<|fe|<n ^ 2 l a fe| 2 } an d the normalization constant, 

Z„ := / ■•• / e" 27r2 ^ fc2|afc|2 daida_i---da„da_ n . 
Jc •/: 

We note that this is a petit, or mixed, canonical ensemble: canonical in H and micro- 
canonical in the mass, i.e., the L 2 norm is constrained to be no more than a constant 
B. The limit of p n as n — > oo is p, which turns out to be the image measure of the 
random Fourier series map: 



ui i — y 



k^O 



2nk 



The gk are independent identically distributed complex Gaussian variables, and the 
image is almost surely in H s for any s < |. And the limit of p n as n — > oo is /x, the 
following weighted Wiener measure: 



(4.4) 



dp = Z 



INUi 



L {||u|ll^s} da o ® dp- 

Here dp can be represented either through the random Fourier series leaving off the 
zero mode, or the formal Weiner measure (with its own normalization Zq): 



da D ®dp = Z^e-ift |a * u|2(fa ]J du(a 



x£T 



Then the invariance of the measures p n under the finite-dimensional flow S n of the 
spectrally truncated NLS (4.3), combined with the uniform regularity of S n in H s and 
the approximation of the NLS flow by S n , all give the following theorem for the limit 
measure: 



Theorem 4.1 (Bourgain [3]). The measure p in (4.4 1 is invariant under the dynamics 
of the periodic ID focusing cubic (p = 3) and quintic (p — 5) NLS. 

This is also parlayed into a global well-posedness result for the continuous NLS: 

Theorem 4.2 (Bourgain [3]). The Cauchy problem for the ID periodic NLS with p 5 
and with random initial data: 

\ id t u = -d xx u - |k| p_1 u 
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is globally well-posed for almost all ao and to satisfying the mass constraint: 

Another advantage of Gibbs measures has been illustrated in recent work on the ID 
cubic NLS with low-regularity initial data on the circle T, namely that in spite of the 
ill-posedness results for the NLS with initial data of lower regularity than L 2 [17] , the 
bad initial data are quite rare, forming a set of vanishing measure with respect to the 
natural Gibbs measure. In fact, for the cubic NLS on the circle, Colliander and Oh 
were able to prove almost-sure local well-posedness for initial data in £f <7 (T), for every 
a > —1/3, and global well-posedness for data in H S (T), for s > —1/12 [T5]. Oh and 
Sulem proved weak continuity of the cubic NLS solution map from L 2 (T) to the space 
of distributions [48] . Both papers relied on a Wick ordering of the nonlinearity, which 
seems to be an appropriate approach for the ID cubic NLS with rough initial data; by 
contrast Molinet proved that the standard solution map (without Wick ordering) from 
L 2 (T) to the space of distributions is not weakly continuous. This approach has also 
been fruitful for almost sure global well-posedness for the periodic derivative NLS [17] . 

4.2. Application of Gibbs measures to solitons. Jordan, Josserand, and Turking- 
ton have established the connection between Gibbs measures and solitons for the ID 
NLS with non-traditional (i.e., not polynomial) nonlinearities [221 ISO]- Consider the 
following NLS on a finite interval [0, L] with periodic boundary conditions: 

(4.5) id t u = -d xx u - g(\u\ 2 )u, 

with focusing and bounded nonlinearity, i.e., <?(0) =0, Va > 0,g'(a) > and g(a) ^ C, 
and as a — > oo,ag'(a) ^ C"; for example a nonlinearity that approximates behavior of 
large wave amplitudes: 

n i2\ M 2 

sM = ' i2 - 

As usual the two conserved quantities are energy: 

H(u) -=\J (\d x it\ 2 - G{\u\ 2 ))dx, 
where G is the primitive of g: G(a) = g(s)ds, and mass: 

M{u) :=\ f \u\ 2 dx. 
2 Jo 

In this ID setting, the energy appearing in the Gibbs measure contains the kinetic 
term which captures the fine-grained fluctuations and the potential term which essen- 
tially ignores those fluctuations. 

Using the same kind of spectral truncation as Bourgain, the mixed ensemble measure 
is formed with a regularization with respect to mass as well as the usual conditioning 
on constrained mass: 

Pf (du n ) := (Z(j3 n , B))- 1 exp (0 [#(«") + v\\u n \\\])l { \\ u ^ B} J] du n . 

They fix the mass B and scale the inverse temperature as (3 n ~ Cn for large n so 
that the mean energy has a finite limit, then rigorously analyze the limit using large 
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deviations theory, comparing the result with numerical simulations |29j . The large 
deviations principle is that for appropriate sets S C L 2 : 

lim — log Pg {u n E B} = - inf I B (u), 
where I B is the large deviations rate function defined for u £ H^: 

I B (U) ~I H ^- E W> -2\\ U W 2 2 = B > 

1 oo, otherwise. 

Here E is the coherent energy function E(N) := mm{H(u) : M(u) = B,u E Hq}. 

What this means is that with high probability the random fields (spectrally truncated 
wavefunctions) u n converge exactly to the ground state solitons of the continuous NLS. 
This is because the minimizers of the large deviations rate function are the states that 
are most probable macroscopically, and all other states have exponentially vanishing 
probability. And since the large deviations rate function is equal to the Hamiltonian 
energy up to a constant, the ground state solitons that solve the traditional variational 
problem are identical with the most probable macroscopic states. 

The ground states also turn out to be stable both macroscopically and microscopi- 
cally: the large deviations principle says that the macrostate stays close to the ground 
states with high probability, even if some of the energy goes into turbulent fluctuations; 
and the typical microstate as sampled from the Gibbs measure P B n is close to a ground 
state in the .ff^-norm. 

Two key things that Jordan, Josserand, and Turkington need in order for their 
approach to work are non-integrability and the absence of wave collapse, so in particular, 
it doesn't work for the cubic NLS in one, two, or three dimensions. It turns out that the 
ID and 3D cases can be handled in the right frameworks, and possibly two dimensions 
as well. 

4.3. Conjecture about a phase transition to soliton-like behavior. Numerical 
simulations led Lebowitz, Rose, and Speer to make a second conjecture about a possible 
phase transition in the invariant measures for the ID NLS: at high temperature and 
low mass they saw uniformly small wavefunctions, but at low temperature and high 
mass there emerged sharply concentrated structures like solitons. 

By contrast, simulations in [10 suggested that there is no phase transition, and 
then Rider, following on the works of McKean and Vaninsky [41j [42j [43] , confirmed 
this for the ID infinite-volume focusing NLS by proving that the thermodynamic limit 
is trivial [52, 53 . Bourgain also studied invariant measures of the ID infinite- volume 
defocusing NLS [S] and of the 2D defocusing NLS [5] (see also the review article [7]), 
as did Tzvetkov [59] , 

These statistical mechanics methods that work so well for the subcritical NLS break 
down at criticality. For instance, Brydges and Slade studied the 2D focusing NLS 
[5] and saw that the natural Gibbs measure construction cannot produce an invariant 
measure for large coupling coefficients of the (critical) nonlinearity. Similarly for 3D 
focusing NLS, the natural construction is not normalizable, and it is thought to be 
impossible to make any reasonable Gibbs construction. 

This is the motivation for the alternative approach explained below, using a spatial 
discretization to construct the finite-dimensional approximation corresponding to the 
focusing cubic NLS in three dimensions. There are still unresolved issues about what 
information can be transferred to the continuum limit, but this approach has turned 
out to be very fruitful, especially as far as how much can be said about properties of 
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the Gibbs measures-as it turns out, a lot more than the previous spectral approaches, 
and even about a phase transition to solitons. 

5. A PHASE TRANSITION TO SOLITONS FOR 3D DISCRETE NLS 

Because constructing Gibbs measures for the 3D NLS appears to be impossible, 
we consider instead the NLS on a spatial discretization of the 3D unit torus [0, l] 3 , 
represented by: 

V = {0,1/L,...,(i-1)/L} 3 , 

so that each vertex of V has degree d — 3, V has n := \V\ — L 3 sites, and the inter-site 
distance is h = 1/L. 

(We consider this interesting special case here for simplicity, but in [16j the results 
actually handle general graphs G — (V, E) that are finite and undirected without self- 
loops, and that are high- dimensional in the sense that nh 2 goes to oo as h — > 0, where 
h > is the distance between two neighboring vertices in G and n is the number of 
sites, n — > oo. This condition holds for the discrete 3D torus defined above, but does 
not hold for the discrete 2D torus in the usual continuum limit.) 

The focusing NLS on the discrete 3D torus is: 

(5.1) ij t f k = -A/ fc - \f k \ 2 f k , 

with discrete Laplacian: 

(5-2) Af k = p J2 Ui - A)- 

j:j~k 

Solutions of (5.1 ) are known to exist for all time [60 , and the dynamics of this system 
conserves the power N(f) := J2kev l/fc| 2 ' or equivalently mass M(f) := n~ 1 N(f), and 
the energy, appropriately normalized by the number of sites n: 

2 



"(f) - ; E 

j,kev 



fk - fj 



n 

k£V 



Then by the Liouville theorem, the natural Gibbs measure d/i := e P H W> Ylkev dfk is 



invariant under the dynamics of the discrete NLS (5.1 ) for any real /3 (inverse temper- 
ature). However, this measure has infinite mass if /3 > 0, so we construct the Gibbs 
measures with a mass cut-off to restrict the system to an allowed mass B: 

(5.3) dfo, B := Z- l e-^ f h {N(mBn} J] df k . 

kev 

Here B is arbitrary and positive, and the normalizing constant, or partition function, 
is Z = Z(f3,B). This cut-off Gibbs measure is still invariant under the NLS dynamics 
because mass is a conserved quantity, and it turns out to be tractable in the limit as 
the grid size h goes to zero, at least in dimensions three and higher. 

Let ^ be a random element of with probability distribution fl: that is, i/> is a 
random wavefunction on V such that for each A C C , 

F(^eA) = Z- 1 [ e-P H Vh {N(f) < Bn} df, 

J A 

where df = JT df k denotes the Lebesgue differential element on C v . In order to 
understand the behavior of the random map ip, we first study the partition function 
Z. The first theorem below shows that if we have a sequence of graphs with n and nh 2 
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P 

e' c =2.455... *" 

Figure 1. The free energy is constant for small inverse temperature 
and starts increasing at the critical threshold. Here mass is normalized, 
B = 1. 



both tending to infinity, the limit of n 1 log Z can be exactly computed for any positive 
/3 and B. 

5.1. The partition function and free energy. The first result is that for positive 
temperature, and for the discrete torus in dimensions three and higher, there exists a 
function F that we construct, such that as n — > oo: 

Z{p,B)~e nF ^> B \ 

More precisely, let g : [2, oo) — » R be the function 



(5.4) ^,».| + ^ 1 _| + v,(l-^l- 

It can be checked that g has a unique real zero that we call 9 C , because g is strictly 
increasing in [2, oo), g(2) < and g(3) > 0. Numerically, 8 C ss 2.455407. Let 



(5.5) F(fi,B):= 



log(BTre) if/3S 2 <6» c , 
log(STre) + g{PB 2 ) if (3B 2 > 6 C . 

(Figure [T] shows a graph of F versus /? when B = 1.) The theorem also gives an explicit 
rate of convergence. 

Theorem 5.1 (Chatterjee-Kirkpatrick [TB]). Suppose /3 > 0. Let e £ (0,1/5). There 
exists a positive constant C depending only on e, f3, B , h and d such that if n > C , 
then 

— - F{B, B) > -Cn- 1 / 5 - C{nh 2 )- 1 
n 

and 

logZ _ < iCn- 1 '^ + Cn~**/* if PB 2 < 6 C , 

n [P ' ) ~ yCn- 1 ' 5 ^ if/3B 2 > 9 C . 

The proof of these asymptotics for the partition function is at the heart of the work 
[16] : large deviations techniques for discrete random variables are used, first obtaining 
upper and lower bounds for lo ^ z and then optimizing which gives the free energy F. 
Then, because F is not differentiable at the critical threshold 9 C , that means there is 
a first-order (discontinuous) phase transition at 9 C . The random map ip behaves quite 
differently in the two phases, and there is an important quantity which is discontinuous 
at the threshold. 



14 



KAY KIRKPATRICK 



5.2. The phase transition and its consequences. To roughly describe this phase 
transition (depicted in figure [2]), when f3B 2 > 9 C there is a single site k € V which 
bears an abnormally large fraction of the total mass of the random wavefunction ip. 
This fraction is nearly deterministic, given by the ratio a/B, where 



(5.6) a = a(J3,B) := ~ + 



2 V PB 2 ' 




Figure 2. Impressionistic pictures of typical wavefunctions in the two 
phases. Left: subcritical noncondensed phase, f3B 2 < 9 C , where the 
heaviest site bears a negligible proportion of the mass compared to 
the total mass N(ip) ~ Bn. Right: supercritical condensed phase, 
/3B 2 > 6 C , where a significant proportion a of the mass concentrates 
at a single site. 



More precisely, if M\(ip) and Ma(V') are the largest and second largest components 
of the vector (\il>k\ 2 )kev> then with high probability Mx(ip) ~ an and Ma^) = o(n). 
Moreover, N(ip) w Bn with high probability. A consequence is that the largest com- 
ponent carries more than half of the total mass: 

IV'fcl 2 a 1 

On the other hand, when /3B 2 < 9 C , then Mi(ip) = o{n), but still N(ip) w Bn. 
Consequently 

l^fcl 2 n 

The formula for a (see also the graph of a in figure [3] as a function of /3 for B = 1) 
shows that a does not tend to zero as (3B 2 approaches 9 C from above, making it clear 
that this is a discontinuous (in physics parlance, first-order) phase transition. 

Also, in the regime f3B 2 < 8 C , the energy density is close to zero. By contrast, when 
j3B 2 > 9 C , the energy density H(ip)/n is strictly negative and approximately equals 
—a 2 , the exceptional site bearing nearly all of the energy of the system. This is because 
the total energy H(ip) is approximately — a 2 n, while the energy at k is, summing over 
just the neighbors j of k: 

1 V-,, ,,2 I/-' 4 



nh 2 



= —a 2 n + o(n), 

the equality by Theorem |5.2[ which makes precise this phase transition: 
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Figure 3. The fraction of mass at the heaviest site jumps from 
roughly zero for small inverse temperature, to roughly .71 at the crit- 
ical threshold. (Here B — 1.) 



Theorem 5.2 (Chatterjee-Kirkpatrick [16] )■ Let V be the discrete torus in dimension 
d ^ 3, let a — a(j3,B) be defined as in ( |5.6| , and let Mi(ip) and M2(ip) be the largest 
and second largest components of {\4>k\ 2 )kev ■ 

(1) First, suppose (3B 2 > 9 C . For any q £ (|, 1), there is a constant C depending 
only on fj, B, d, and q such that if n > C, then with probability > 1 — e~ n l c , 



(5.7) 



n 
n 



< Cr 



-(i-«)/4 



< Cn- (1 - q ^ 4 , and 



-B <Cn-^' 2 



(2) Next, suppose j3B 2 < 9 C . For any q G (j|,l), there is a constant C depending 
only on (3, B, d, and q such that whenever n > C, with probability > \ ~e~ nq l c 



(5.8) 



iV(V>) 



B 



< n 



-(1-9) 



and 



Mi WO 



(3) Finally, if (3B 2 — 8 C and q £ (y|, 1), then there is a constant C depending only 



on (3, B, d, p and q such that whenever n > C , with probability > 1 — e 
either (5.7) or (5.8 1 holds. 



q /C 



The proof proceeds by a multi-step approximation: first the Gibbs measures are 
approximated by uniform measures on appropriate spaces (intersections of I 2 and I 
spheres), and then the uniform measures are approximated by i.i.d. complex Gaussian 
random variables on Z 4 annuli. Then the i.i.d. Gaussians must concentrate on the parts 
of the Z 4 annuli closest to the origin, i.e., where one coordinate is large and the others 
are uniformly small. 

One consequence of this theorem is confirmation of the conjecture of Flach, Kladko, 
and MacKay, because the energies of the localized states are bounded away from zero, 
made clear in figure [4] where energy is plotted as a function of f3. 
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Figure 4. Localized states occur only when the energy is nonzero, 
and in fact bounded away from zero, as in the physics conjecture. 



Another straightforward consequence of this theorem is that typical discrete wave- 
functions above the critical threshold blow up in the discrete H 1 norm: 

2 



1 



kev 



\fk 



E 

(kj)eE 



fk ~ fj 



It also turns out that the discrete H 1 norm also diverges even for (3B 2 < 6 C , just with 
different rates: 



Theorem 5.3 (Chatterjee-Kirkpatrick [H]). Suppose the context of Theorem 5.2 holds, 
in dimension d ^ 3, with n > C. Let ip = (ipk)kev be a discrete wavefunction picked 
randomly from the invariant probability measure ft defined in (5.3 1. If f3B 2 > 9 C , then 
there is a positive constant c depending only on j3, B, d and p such that PfllV'llrfi < 
cn p ) < e~ n whenever n > 1/c. On the other hand, if f3B 2 < 6 C , then similarly: 
Si < cVdn p ) < e" 



-dn a 



Since the measure ft of (5.3) is invariant for the discrete NLS (5.1), one may expect 



that if the initial data comes from ft, localized modes will continue to exist for all time. 
The question is whether the mode jumps around or stays in one place (in which case 
we have a standing or stationary wave, also called a discrete breather). The following 
theorem uses the statistical equilibrium result to deduce a dynamical result: the same 
site continues to be the localized mode for an exponentially long period of time. 



Theorem 5.4 (Chatterjee-Kirkpatrick [16]). 



j3B > 9 C and the function a as defined in (5.6) 



Suppose the previous context holds with 
Let ip(t) 



be a discrete 



(ipk(t))kev 

wavefunction evolving according to (5.1), where the initial data ij}(0) is picked randomly 

'4 
*5 ' 



from the invariant probability measure ft defined in (5.3). Then for any q € (|,1), 
there is a constant C depending only on f3, B, d, and q such that if n > C ', then with 
probability > 1 — e~ n? '/ c the inequalities (5.7) hold for if)(t) for all < t < e™ 9 /* 7 



and moreover there is a single k e V such that the maximum of \ipj(t)\ * s attained at 
y — x for all Q < t < e nq ^ c . In particular, ip(t) is approximately a standing wave with 
localized mode at k for an exponentially long time. 

This dynamical stability result combines PDE and probabilistic methods. First, with 
high probability the wavefunction is condensed at time t — 0, and then the invariance 
of the Gibbs measure implies that the wavefunction is condensed for an overwhelming 
proportion of times up to the exponential time in the theorem. If the soliton were to 
jump from one site to another, then it would have to happen in a tiny window of time; 



but that would result in a wavefunction that violates basic NLS estimates, so in fact 
the soliton must stay put. 



In other words, it is not only possible, but typical, for solutions of ( 5.1 ) to have unique 
stable localized modes for exponentially long times if the initial inverse temperature or 
mass are above a threshold. Earlier results for localized modes were only for existence, 
and only for large mass, i.e., tending to infinity (e.g., [60J, while Theorem 5.4 proves 
typicality for finite mass. 

Our final result is about the probability distribution of the individual coordinates of 
a random map if) picked from the measure fip,B- We can give a precise description of the 
distribution for small collections of coordinates: In both phases, the joint distributions 
are approximately complex Gaussian vectors with the appropriate variance (B below 



the threshold and B — a above it, with a defined in (5.6)), provided the collection is 
sufficiently small compared to the number of sites in the grid and provided a symmetry 
assumption is satisfied. 

The symmetry assumption on G is that G is translatable by some group S: where E 
is a group of automorphisms of G such that: 1) |S| = n, and 2) no element of £ except 
the identity has any fixed point. For example, the discrete torus is translatable by the 
group of translations. 

Theorem 5.5 ( Chatter jee-Kirkpatrick 16 ). Suppose the graph G is translatable by 
some group of automorphisms, d 3, and ijj is a random wavefunction picked according 
to the measure fi. Take any m distinct points X\, . . . , x m G V . Let cf> = (<f>i, . . . , </>„,.) 
be a vector of i.i.d. standard complex Gaussian random variables. If f3B 2 < 8 C , then 
there is a constant C > depending only on /?, B, and d such that if n > C , then for 
all Borel sets U C C"\ 

\¥{B- l ' 2 (^ Xl , . . . , % n ) G U) - P(0 G U)\ < mrr l ' c . 

If f3B 2 > 9 C , the result holds after B^ 1 / 2 is replaced with (B — a) -1 / 2 where a = a(/3, B) 
is defined in (5.6), and the error bound is changed to vr?n~ x l c . 



This says that below the critical threshold, the wavefunctions look like white noise, 
but above the threshold they look like a single large coordinate and white noise else- 
where. 

As for the future, one direction to push these kinds of results is into the two- 
dimensional setting. There the nonlinear and kinetic energies arc in a delicate balance, 
and the wavefunction can't be approximated by i.i.d. Gaussians; instead one must 
consider correlations between sites, producing what is called a discrete Gaussian free 
field. Another direction is to study the Gibbs measures for the fractional NLS, which 
will have non-local behavior coming from the underlying Levy processes. There is also 
a possibility that this statistical mechanics approach could be useful for the soliton res- 
olution conjecture: that in certain regimes, solutions generically tend asymptotically 
towards a soliton plus the small fluctuations of a radiation component. 
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